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Abstract 

The time irreversibility problem is the dichotomy of the reversible microscopic dy- 
namics and the irreversible macroscopic physics. This problem was considered by Boltz- 
mann, Poincare, Bogolyubov and many other authors and though some researchers 
claim that the problem is solved, it deserves a further study. In this paper an attempt 
is performed of the following solution of the irreversibility problem: a formulation 
of microscopic dynamics is suggested which is irreversible in time. In this way the 
contradiction between the reversibility of microscopic dynamics and irreversibility of 
macroscopic dynamics is avoided since both dynamics in the proposed approach are 
irreversible. 

A widely used notion of microscopic state of the system at a given moment of time 
as a point in the phase space and also a notion of trajectory and microscopic equation 
of motion does not have an immediate physical meaning since arbitrary real numbers 
are non observable. In the approach presented in this paper the physical meaning is 
attributed not to an individual trajectory but only to a bunch of trajectories or to the 
distribution function on the phase space. 

The fundamental equation of the microscopic dynamics in the proposed "functional" 
approach is not the Newton equation but the Liouville equation for the distribution 
function of the single particle. Solutions of the Liouville equation have the property of 
derealization which accounts for irreversibility. It is shown that the Newton equation 
in this approach appears as an approximate equation describing the dynamics of the 
average values of the position and momenta for not too long time intervals. Corrections 
to the Newton equation are computed. 
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1 Introduction 



The time irreversibility problem is the problem of how to explain the irreversible behaviour 
of macroscopic systems from the time-symmetric microscopic laws. The problem has been 
discussed by Boltzmann, Poincare, Bogolyubov, Kolmogorov, von Neumann, Landau, Pri- 
gogine, Feynman and many other authors [1] - [19] and it deserves a further study. 

In particular, in works by Poincare [2], Landau and Lifshiz [5], Prigogine [TO] . Ginzburg 
[14] . Feynman [16] it is stressed that the irreversibility problem is still an open problem. 
Poincare [2] said that perhaps we will never solve the irreversibility problem. Landau and 
Lifshiz write about the principle of increasing entropy [5]: "Currently it is not clear whether 
the law of increasing entropy can be in principle derived from classical mechanics." Lan- 
dau speculated that to explain the second law of thermodynamics one has to use quantum 
mechanical measurement arguments. 

From the other side Lebowitz [17], Goldstein [18] and Bricmont [19] state that the irre- 
versibility problem was basically solved already by Boltzmann by using his notion of macro- 
scopic entropy and the probabilistic approach. 

The microscopic mechanical description of a system assumes that the state of the system 
at a given moment of time is represented by a point in the phase space with an invariant 
measure and the dynamics of the system is described by a trajectory in the phase space, see 
[5J [201 ED E2J E31 EI] • It is assumed that the microscopic laws of motion are known (Newton 
or Schrodinger equations) and there is a problem of derivation from them the macroscopic 
(Boltzmann, Navier-Stokes,...) equations, see for example, [3 [8]. 

There are well known critical remarks by Loschmidt and Poincare and Zermelo on the 
Boltzmann approach to the irreversibility problem and the .ff-theorem. Loschmidt remarked 
that from the symmetry of the Newton equations upon the reverse of time it follows that to 
every motion of the system on the trajectory towards the equilibrium state one can put into 
correspondence the motion out of the equilibrium state if we reverse the velocities at some 
time moment. Such a motion is in contradiction with the tendency of the system to go to 
the equilibrium state and with the law of increasing of entropy. 

Then, there is the Poincare recurrence theorem which says that a trajectory of a bounded 
isolated mechanical system will be many times come to a very small neighborhood of an initial 
point. This is also in contradiction with the motion to the equilibrium state. This is the 
Poincare-Zermelo paradox. 

Boltzmann [25] gave the following answer to the Loschmidt argument: "We do not have 
to assume a special type of initial condition in order to give a mechanical proof of the 
second law, if we are willing to accept a statistical viewpoint. While any individual non- 
uniform state (corresponding to low entropy) has the same probability as any individual 
uniform state (corresponding to high entropy), there are many more uniform states than 
non-uniform states. Consequently, if the initial state is chosen at random, the system is 
almost certain to evolve into a uniform state, and entropy is almost certain to increase." 

So, the answer by Boltzmann to the objection of Loschmidt was that, firstly, the proba- 
bilistic considerations has been involved, and secondly, he argued that with the overwhelm- 
ing probability the evolution of the system will be occur in the direction of flow of time, 
corresponding to the increasing entropy, since there are many more uniform states than non- 
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uniform states. The answer by Boltzmann to the Poincare - Zermelo objection was in the 
pointing out the extremely long Poincare recurrence time. 

These Boltzmann's responses are not very convincing, from our point of view, despite 
their vigorous support in recent works [TT], [THJ [TH] - Involvement of probability considerations 
alone does not clarify the issue of irreversibility, because if there is symmetry in relation to 
the direction of time, it remains unclear why the evolution in one direction is more likely 
than the other. 

Then, the argument, that there are many more uniform states than non- uniform states 
does not clarify the issue of the dynamical evolution since the dynamics does depend on 
the form of the potential energy between particles and for many potentials the argument is 
simply wrong. Therefore this general Boltzmann's argument does not give a real insight to 
the irreversibility problem. 

Actually, Boltzmann in [25] considered "a large but not infinite number of absolutely 
elastic spheres, which move in a closed container whose walls are completely rigid and likewise 
absolutely elastic. No external forces act on our spheres." Even for this simple model it is 
very difficult to make the Boltzmann argument convincing, i.e. to get a mathematical result, 
see j2HI2g. 

Further, an indication to the extremely long Poincare recurrence time does not remove the 
contradiction between microscopic reversibility and macroscopic irreversibility, and moreover 
no clear mechanism for relaxation to equilibrium is presented. 

Lebowitz advanced p2], following to Boltzmann, the following arguments to explain 
irreversibility: a) the great disparity between microscopic and macroscopic scales, b) a low 
entropy state of the early universe, and c) the fact that what we observe is the behaviour of 
systems coming from such an initial state - not all possible systems. 

From our viewpoint these arguments do not lead to explanation of irreversibility even 
though it is said in [T7] that "common alternative explanations, such as those based on the 
ergodic or mixing properties of probability distribution... are either unnecessary, misguided 
or misleading" . 

Boltzmann proposed that we and our observed low-entropy world are a random fluctu- 
ation in a higher-entropy universe. These cosmological considerations of the early universe 
might be entertaining but they should be related with the modern Friedmann [21^ |2"T] grav- 
itational picture of the Big Bang and, what is most important, there is no evidence that 
the irreversible behaviour of gas in a box is related somehow with conditions in the early 
universe 14 billions years ago. 

Notice that in [28j it is shown that the Hawking black hole information paradox is a 
special case of the irreversibility problem. 

Goldstein said in [T5]:"The most famous criticisms of Boltzmann's later works on the 
subject have little merit. Most twentieth century innovations - such as the identification 
of the state of a physical system with a probability distribution p on its phase space, of its 
thermodynamic entropy with the Gibbs entropy of p, and the invocation of the notions of 
ergodicity and mixing for the justification of the foundations of statistical mechanics - are 
thoroughly misguided." 

And then: "This use of ergodicity is thoroughly misguided. Boltzmann's key insight was 
that, given the energy of a system, the overwhelming majority of its phase points on the 
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corresponding energy surface are equilibrium points, all of which look macroscopically more 
or less the same." 

The Boltzmann argument about "the overwhelming majority" (i.e. "many more uniform 
states") was discussed above. Moreover, the main point of the current paper is that we shall 
use the probability distribution and the Liouville equation not only in statistical mechanics 
but also in classical mechanics, even for a single particle in empty space. 

A powerful method for obtaining kinetic equations from the Newton - Liouville equations 
was developed by Bogolyubov [3]. He has considered infinite number of particles in infinite 
volume and postulated the condition of weakening of initial correlations between particles 
in the distant past, through which the irreversibility entered into the equation for the dis- 
tribution functions, as well as using a formal expansion in powers of density, which leads to 
divergences. 

Poincare considered the model of free motion of gas particles in a box with reflecting walls 
and showed that for solutions of the Liouville equation in this model there is, in some sense, 
an irreversible diffusion [2S]. This result of Poincare was introduced to modern scientific 
literature by Kozlov, see [IT], where the result of Poincare was significantly strengthened 
and consolidated. In the works of Kozlov a method of the weak limit in the nonequilibrium 
statistical mechanics has been developed, and, in particular, it was proved that for some 
models the system in the sense of weak convergence tends to one and the same limit in 
the past and in the future [TT1 [T2] . The method of the weak limit of [TT1 [T2] had a signif- 
icant influence to the formulation of the approach to the problem of irreversibility through 
functional formulation of classical mechanics. 

Note that the stochastic limit [13] gives a systematic method for investigation of irre- 
versible processes. 

Questions about the increase of the fine and coarse entropies are discussed in [301 El El 

LTl El EH LT23- 

In this paper we attempt to suggest the following approach to the irreversibility prob- 
lem and to paradoxes of Loschmidt and Poincare - Zermelo: we propose a formulation of 
microscopic dynamics which is irreversible in time. Thus the contradiction between micro- 
scopic reversibility and macroscopic irreversibility of the dynamics disappears, since both 
microscopic and macroscopic dynamics in the proposed approach are irreversible. 

Note that the conventional widely used concept of the microscopic state of the system 
at some moment in time as the point in phase space, as well as the notion of trajectory 
and the microscopic equations of motion have no direct physical meaning, since arbitrary 
real numbers not observable (observable physical quantities are only presented by rational 
numbers, cf. the discussion of concepts of space and time in [331 EH ESI ESI EZl EBJ ESI BE])- 

In the proposed "functional" approach, the physical meaning is attached not to a single 
trajectory, but only to a "beam" of trajectories, or the distribution function on phase space. 
Individual trajectories are not observable, they could be considered as "hidden variables", if 
one uses the quantum mechanical notions, see [HI 02]. 

The fundamental equation of the microscopic dynamics of the proposed functional proba- 
bilistic approach is not Newton's equation, but a Liouville equation for distribution function. 
It is well known that the Liouville equation is used in statistical mechanics for the descrip- 
tion of the motions of gas. Let us stress that we shall use the Liouville equation for the 
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description of a single particle in the empty space. 

Although the Liouville equation is symmetric in relation to the reversion of time, but 
his solutions have the property of delocalization, that, generally speaking, can be interpreted 
as a manifestation of irreversibility. It is understood that if at some moment in time the 
distribution function describes a particle, localized to a certain extent, then over time the 
degree of localization decreases, there is the spreading of distribution function. Delocalization 
takes place even for a free particle in infinite space, where there is no ergodicity and mixing. 

In the functional approach to classical mechanics we do not derive the statistical or 
chaotic properties of deterministic dynamics, but we suggest that the Laplace's determinism 
at the fundamental level is absent not only in quantum, but also in classical mechanics. 

We show that Newton's equation in the proposed approach appears as an approximate 
equation describing the dynamics of the average values of coordinates and momenta for not 
too long time. We calculate corrections to Newton's equation. 

In the next section the fundamentals of the functional formulation of classical mechanics 
are presented. Sections 3 and 4 deal with the free movement of particles and Newton's 
equation for the average coordinates. Comparison with quantum mechanics is discussed in 
Section 5. General comments on the Liouville and Newton equations are contained in section 
6. Corrections to the Newton equation for a nonlinear system are calculated in Section 7. 
Reversibility of motion in classical mechanics and irreversibility in the functional approach 
to the mechanics discussed in section 8. The dynamics of the classical and quantum particle 
in a box and their interrelationships are summarized in section 9. 

2 States and Observables in 

Functional Classical Mechanics 

Usually in classical mechanics the motion of a point body is described by the trajectory in 
the phase space, i.e. the values of the coordinates and momenta as functions of time, which 
are solutions of the equations of Newton or Hamilton. 

Note, however, that this mathematical model is an idealization of the physical process, 
rather far separated from reality. The physical body always has the spatial dimensions, 
so a mathematical point gives only an approximate description of the physical body. The 
mathematical notion of a trajectory does not have direct physical meaning, since it uses 
arbitrary real numbers, i.e. infinite decimal expansions, while the observation is only possible, 
in the best case, of rational numbers, and even them only with some error. Therefore, in the 
proposed "functional" approach to classical mechanics, we are not starting from Newton's 
equation, but with the Liouville equation. 

Consider the motion of a classical particle along a straight line in the potential field. The 
general case of many particles in the 3-dimensional space is discussed below. Let (q,p) be 
co-ordinates on the plane M 2 (phase space), t G K is time. The state of a classical particle 
at time t will be described by the function p = p(q,p,t), it is the density of the probability 
that the particle at time t has the coordinate q and momentum p. 

Note that the description of a mechanical system with the help of probability distribution 
function p = p(q,p, t) does not necessarily mean that we are dealing with a set of identically 
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prepared ensemble of particles. Usually in probability theory one considers an ensemble of 
events and a sample space [551 ESI [57] . But we can use the description with the function 
P — p{QjP,t) also for individual bodies, such as planets in astronomy (the phase space in 
this case the 6-dimensional). In this case one can think on the "ensemble" of different 
astronomers which observe the planet, or on the "ensemble" of different models of behaviour 
of a given object. Actually, it is implicitly always dealt with the function p = p(q,p, t) which 
takes into account the inherent uncertainty in the coordinates and momentum of the body. 
An application of these remarks to quantum mechanics will be discussed in a separate work. 

The specific type of function p depends on the method of preparation of the state of a 
classical particle at the initial time and the type of potential field. When p = p(q,p,t) has 
sharp peaks at q = qo and p = Po, we say that the particle has the approximate values of 
coordinate and momentum q and p . 

Emphasize that the exact derivation of the coordinate and momentum can not be done, 
not only in quantum mechanics, where there is the Heisenberg uncertainty relation, but also 
in classical mechanics. Always there are some errors in setting the coordinates and momenta. 
The concept of arbitrary real numbers, given by the infinite decimal series, is a mathematical 
idealization, such numbers can not be measured in the experiment. 

Therefore, in the functional approach to classical mechanics the concept of precise trajec- 
tory of a particle is absent, the fundamental concept is a distribution function p = p(q,p,t) 
and 5-function as a distribution function is not allowed. 

We assume that the continuously differentiable and integrable function p = p(q,p,t) 
satisfies the conditions: 



The formulation of classical mechanics in the language of states and observables is considered 
in [131FHIB5] . The functional approach to classical mechanics differs in the following respects. 
Because the exact trajectory of a particle in the functional approach does not exist, then the 
function p = p(q,p,t) can not be an arbitrary generalized function, it is the usual function 
of class L 1 (R 2 ), or even continuously differentiable and integrable function. 

In addition, the motion of particles in the functional approach is not described directly 
by the Newton (Hamilton) equation. Newton's equation in the functional approach is an ap- 
proximate equation for the average coordinates of the particles, and for non-linear dynamics 
there are corrections to the Newton equations. 

As is known, the mathematical description of a moving fluid or gas is given by means of 
the density distribution functions p(q,t), as well as the velocity v(q,t) and pressure p(q,t), 
see, for example, [16]. Let the function p(q,p,t) describes a particle, as proposed in the 
functional formulation of classical mechanics, and we set p c (q,t) = f p(q,P,t)dp. We could 
ask the question can we determine by the form of functions p(q,t) and p c (q,t) whether we 
are dealing with a continuous medium or with a particle? The general answer is the follow- 
ing: functions p(q,t) and p c (q,t) satisfy different equations (the Navier-Stokes or Liouville 
equation) and different conditions of normalization. 

Note, however, that if an error in determining the coordinates and momentum of particles 
is large enough, it really is not so easy to determine, we have a case of, say, a fast-moving 
particle in a box with reflecting walls, either a gas of particles. 
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If / = f(q,p) is a function on phase space, the average value of / at time t is given by 
the integral 

/(*)= / f(q,p)p(q,P,t)dqdp. (2) 



In a sense we are dealing with a random process £(£) with values in the phase space. Motion 
of a point body along a straight line in the potential field will be described by the equation 

dp _ P dp | dV(g) dp 
dt m dq dq dp 

Here V(q) is the potential field and mass m > 0. 

Equation ([3]) looks like the Liouville equation which is used in statistical physics to 
describe a gas of particles but here we use it to describe a single particle. 

If the distribution po(q,p) for t = is known, we can consider the Cauchy problem for 
the equation (J3]): 

p\t=o = po(q,p) ■ (4) 
Let us discuss the case when the initial distribution has the Gaussian form: 

1 (q— go) 2 (p-po) 2 

Po{q,P) = — r e a ' 2 e " 2 ■ ( 5 ) 
nab 

At sufficiently small values of the parameters a > and b > the particle has coordinate and 
momentum close to the go and po- For this distribution the average value of the coordinates 
and momentum are: 

q = J qpo(q,p)dqdp = q , p = J pp (q,p)dqdp = p , (6) 

and dispersion 

Aq 2 = {a~^W=^a 2 , Ap 2 =J^W=h 2 . (7) 

3 Free Motion 

Consider first the case of the free motion of the particle when V = 0. In this case the 
equation ([3]) has the form 

dp = _p_dp 
dt m dq 
and the solution of the Cauchy problem is 

P 

p{q,p,t) = po{q — t,p). (9) 

m 

Using expressions (jSJ), ©, 
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we get the time dependent distribution of coordinates: 

/I _ g Q — £2.£) 2 
p{q,p,t)dp= exp{ ™ }, (11) 

while the distribution of momenta is 

/\ (p-pq) 2 
p(q,p,t)dq = -j=^e ~~ P 1- . (12) 

Thus, for the free particle the distribution of the particle momentum with the passage of 
time does not change, and the distribution of the coordinates change. There is, as one says 
in quantum mechanics, the spreading of the wave packet. From ( TTTj) it follows that the 
dispersion Aq 2 increases with time: 

A, 2 (t) = i(a 2 + ^). (13) 
2 m z 

Even if the particle was arbitrarily well localized (a 2 is arbitrarily small) at t — 0, then 
at sufficiently large times t the localization of the particle becomes meaningless, there is a 
delocalization of the particle. 



4 Newton's Equation for the Average Coordinate 

In the functional approach to classical mechanics there is no ordinary picture of an individual 
trajectory of a particle. The starting equation is the dynamic equation ([3]) for the distribution 
function, rather than the Newton equation. 

What role can play the Newton equation in the functional approach? We show that 
the average coordinate for the free particle in the functional approach satisfies the Newton 
equation. Indeed, the average coordinate and momentum for the free particles have the form 

<l(t) = J qpc(q, t)dq = q + ^t, p(t) = J pp m (p, t)dp = p . (14) 

Hence we get 

^?(*) = 0. ( 15 ) 

i.e. we have Newton's equation for the average coordinates. 

We also have Hamilton's equations for the average values of the coordinate and momen- 
tum: 

9=-o^> P=— (16) 
op aq 

where the Hamiltonian H = H(q,p) for the free particle has the form 

E=^-. (17) 
2m v ' 
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Note that in the functional mechanics the Newton equation for the average coordinates is 
obtained only for the free particle or for quadratic Hamiltonians with a Gaussian initial 
distribution function. For a more general case there are corrections to Newton's equations, 
as discussed below. 

We discussed the spreading of Gaussian distribution functions. Similar results are ob- 
tained for the distribution functions of other forms, if they describe in some sense localized 
coordinates and momenta at the initial time. 



5 Comparison with Quantum Mechanics 

Compare the evolutions of Gaussian distribution functions in functional classical mechanics 
and in quantum mechanics for the motion of particles along a straight line. The scene of work 
for the functional classical mechanics is L 2 (IR 2 ) (or L X (IR 2 )), and for quantum mechanics - 
L 2 (IR 1 ). 

The Schrodinger equation for a free quantum particle on a line reads: 

dt 2m dx 2 

Here ip = ip(x,t) is the wave function and h is the Planck constant. The density of the 
distribution function for the Gaussian wave function has the form (see, for example [17]) 

p q (x,t) = mx,t)f = 1 ex P {- iX ~ X ° ~ h2 f^ } ■ ( 19 ) 

a 2 + W (a 2 + 

We find that the distribution functions in functional classical and in quantum mechanics 
( ITT]) and f|T9|) coincide, if we set 

a 2 b 2 = h 2 . (20) 

If the condition (|2"U1) is satisfied then the Wigner function W(q,p,t) [IB] for ip corresponds 
to the classical distribution function (fTUj) , W(q,p, t) = p(q,p,t) . 

The problem of spreading of the quantum wave packet in dealing with the potential 
barrier is considered in [49J. 

Gaussian wave functions on the line are coherent or compressed states. The compressed 
states on the interval are considered in 



6 Liouville Equation and the Newton Equation 

In the functional classical mechanics the motion of a particle along the stright line is described 
by the Liouville equation ([3]). A more general Liouville equation on the manifold T with 
coordinates k ) has the form 

| + t^')-0. (21) 

i=l 



9 



Here p = p(x,t) is the density function and v = v(x) = (v 1 , ...,v k ) - vector field on T. The 
solution of the Cauchy problem for the equation fT2~TT) with initial data 

p\t=o = Po{x) (22) 

might be written in the form 

P(x,t) = po(<P-t(x)) . (23) 
Here <ft{x) is a phase flow along the solutions of the characteristic equation 

x = v(x) . (24) 

In particular, if k = 2n, and M = M n is a smooth manifold, the phase space T = T*M is a 
cotangent bundle, H = H(q,p) is a Hamiltonian function on T, then the Liouville equation 
has the form 

The Liouville measure dp = dqdp is invariant under the phase flow iff. 

Classical dynamical system in the functional approach to mechanics is a stochastic pro- 
cess £(t) = £(t; q,p) = <£t{q,p) which takes values in T and with the probabilistic measure 
dP(q,p) = p (q,p)dqdp. Correlation functions have the form 

<6i(*i) •■■&.(*-) >= J £h(ti, q,p)~4i s (t s ; q,p)po(q,p)dqdp . (26) 

Here ix, i s = 1, k. 

It is assumed usually that the energy surfaces {H = const} are compact. 

A system from iV particles in the 3-dimensional space has the phase space 1R 67V with 

coordinates g = (qi, ...,q7v), P = (pi, Pn), Qi = (ql, qf, qf), Pi = (p},Pi,Pi), i=l,-,N 
and it is described by the Liouville equation for the function p = p(q, p, t) 



dp _ y> ( dV[q) dp p? dp 



dp •s^.dV(q) dp p" dp 

dt ^ dqf dpf rrii dq\ 

Here summation goes on % = 1, N, a = 1,2, 3. The characteristics equations for (j2"7j) are 
Hamilton's equations 

dH . dH 

dpi dqf 



^2 



where the Hamiltonian is 



Emphasize here again that the Hamilton equations (1281) in the current functional approach 
to the mechanics do not describe directly the motion of particles, and they are only the 
characteristic equations for the Liouville equation fT271) which has a physical meaning. The 
Liouville equation ff27|) can be written as 



f = {#>P}> (3°) 
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where the Poisson bracket 

Criteria for essential self-adjointness of the Liouville operator in the Hilbert space L 2 (R 6Ar ) 
are given in [51] . 

7 Corrections to Newton's Equations 

In section 4, it was noted that for the free particle in the functional approach to classical 
mechanics the averages coordinates and momenta satisfy the Newton equations. However, 
when there is a nonlinear interaction, then in functional approach corrections to the Newton's 
equations appear. 

Consider the motion of a particle along the line in the functional mechanics. Average 
value / of the function on the phase space / = f(q,p) at time t is given by the integral (j2J) 

f(t) =< f(t) >= J f(q,p)p(q,p,t)dqdp. (32) 

Here p(q,p,t) has the form ( 1231) 

p(q,P,t) = po(ip- t (q,p)) ■ (33) 

By making the replacement of variables, subject to the invariance of the Liouville measure, 
we get 

</(*)>= / f{q,P)p(q,P,t)dqdp= / f(M ( l'P))po(q,p)dqdp. (34) 



Let us take 

Po(q,p) = S £ (q- q )8 e (p-p ) , (35) 

where 

S e (q) = ^e"^ 2 , (36) 

q E R, e > 0. 

Let us show that in the limit e — > we obtain the Newton (Hamilton) equations: 

\im<f(t)>=f(<p t (q ,p )). (37) 

Proposition 1. Let the function f(q,p) in the expression (SB) be continous and inte- 
grate, and po has the form (E3J). Then 



So / f^P)P^Pi l ) d( l d P = f(ft(qo,Po)) ■ (3* 
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Proof. Functions 5 e (q) form a 5-sequence in D (R) [52]. Hence we obtain 

Jim y f((q,p))p(q,p,t)dqdp = Jim ^ f(ip t (q,p))6 e {q - q )5 t (p - p ) = f{ip t {qo,Po)) , (39) 

that was required to prove. 

Now calculate the corrections to the solution of the equation of Newton. In functional 
mechanics consider the equation, see , 

(40) 

Here A is a small parameter and we set the mass m — 1. The characteristic equations have 
the form of the following Hamilton (Newton) equations: 

p(t) + Xq(t) 2 = , q(t) =p(t). (41) 

Solution of these equations with the initial data 

g (0) = q, q(0) = p (42) 

for small t has the form 

(q(t),p(t)) = <p t (q,p) = (q + pt- ^q 2 t 2 + p - XqH + ...) (43) 
Use the asymptotic expansion 8 e (q) in D (R) for e -> 0, compare [HI [53] : 

4(g) = %) + j5"(g) + -., (44) 

then for e — > we obtain corrections to the Newton dynamics: 

/A e 2 
(q + pt- -q 2 t 2 + ...)[% - g ) + -<y"(g - g ) + •••] (45) 

e 2 A A 

•P(P - Po) + J 5 "(P ~ Po) + -)dqdp = q + p t - -q 2 t 2 - -e 2 t 2 . 

Denoting the Newton solution 

<?Newton(£) = <?0 + p t - ^qfo 2 , 

we obtain for small e, t and A: 

< q(t) >= ^NewtonW ~ J^t 2 . (46) 

Here —^e 2 t 2 is the correction to the Newton solution received within the functional 
approach to classical mechanics with the initial Gaussian distribution function. If we choose 
a different initial distribution we get correction of another form. 
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We have proved 

Proposition 2. In the functional approach to mechanics the first correction at e to the 
Newton dynamics for small t and A for equation has the form 

Note that in the functional approach to mechanics instead of the usual Newton equation 



m^q(t) = F(q) , (47) 



where F(q) is a force, we obtain 



d 2 

m ^2 <?(*) >=< F (<l)( t ) > ■ 



Indeed, multiplying the equation 

g = _JL^_ F( ,)jfc (49) 
at m oq op 

by q and making integration over p and q and then integrating by parts, we get 

Similarly, multiplying the equation (TJ9]) by p and integrating on p and g and then integrating 
by parts, we get 

j t < p(t) >=< F(q)(t) > , (51) 

which gives (H5|) . 

The task of calculating the corrections at e for Newton's equation for mean values is 
similar to the problem of calculating semiclassical corrections in quantum mechanics [5] 



8 Time Reversal 

8.1 Reversibility in classical mechanics 

Let us present a famous discourse which proves reversibility of the dynamics in classical 
mechanics. From the symmetry of Newton's equations upon the replacement the time t for 
—t it follows that if in the system there exists some motion, then it is possible also the 
reverse motion, i.e. such motion, in which the system passes same states in the phase space 
in the reverse order. Indeed, let the function x(t) satisfies the Newton equation 

x(t) = F(x{t)) (52) 

with initial data 

x(0) = xq, x(0) = vq. (53) 
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We denote the corresponding solution by 



x(t) = <&(t;x ,v ) . 

We fix T > and let us reverse the motion of the particle at some moment in time T by 
reversing its velocity, i.e. let us consider the solution y(t) of the Newton equation 

y(t) = F(y(t)) (54) 

with the following initial data: 

y(0) = x(T), y(Q) = -x(T). (55) 

Then it is easy to see that at the time moment T we get 

y{T) = x , y{T) = -v , (56) 

i.e. the particle comes back to the initial point with the inverse velocity. To prove the 
relation (1561) it is enough to note that the solution of equation (Ell) with initial data ( 1551) has 
the form 

y (t) =®(T-t; x 0l v ) 

and use the relations ( f53i) . 

Let us notice that these arguments about reversibility of motion in the classical mechanics 
used not only symmetry of the Newton equation concerning time reversibility, but also the 
fact that a state of the particle in the classical mechanics at some instant of time is completely 
characterized by two parameters - co-ordinate x and speed v. Reversibility of the motion in 
classical mechanics means reversibility of the motion along a given trajectory. 

As it was discussed above, the notion of an individual trajectory of a particle has no 
physical sense. In reality we deal with a bunch of trajectories or probability distribution. In 
the functional classical mechanics the state of the particle is characterized not by the two 
numerical parameters, but the distribution function p = p(q,p, t). In the following subsection 
it will be shown, how it leads to derealization and irreversibility. 



8.2 Irreversibility in the functional mechanics 

The considered reversibility of motion in classical mechanics deals with the individual tra- 
jectory. In the functional mechanics the concept of the individual trajectory of the particle 
has no direct physical sense. Instead, the state of the particle is described by the distribution 
function p = p(q,p,t) which satisfies the Liouville equation (j3J) 

dp _ P dp | dV(g) dp 
dt m dq dq dp 

The Liouville equation is invariant under the replacement t to —t: if p = p{q,p,t) is the 
solution of the equation ( |57l) . then a(q,p,t) = p(q, —p, —t) - also its solution. However this 
symmetry does not mean reversibility of the motion of a particle in the functional approach 
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to mechanics, since the state of the particle is described there by the distribution function 
and the phenomenon of derealization takes place. 

In this way we obtain an answer to the arguments of Loschmidt and Poincare - Zermelo. 
Indeed, to reverse the particle motion at the time moment t = T as it is proposed in the 
Loschmidt argument, it is necessary make the co-ordinate and momentum measurement. But 
it will change the distribution p(q,p,T). Further, it is necessary to prepare such condition 
of the particle that its evolution back in time would lead to the initial distribution p that 
is difficult since the derealization takes place. We will need something even better than 
Maxwell's demon. 

For free particle the derealization leads to the increasing of dispersion Aq 2 with time 

A q \t) = 1 -(a 2 + b ^). (58) 
2 m z 

Notice that the similar phenomena takes place for the Brownian motion B(t) which has 
dispersion t [551 GS] ■ 

Concerning the Zermelo argument related with the Poincare recurrence theorem we note 
that this argument can not be applied to the functional mechanics because this argument 
is based on the notion of individual trajectory. In the functional mechanics the state of 
the system is characterized by the distribution function and here the mean values might 
irreversibly tend to some limits without contradiction with the Poincare theorem as it will 
be shown in the next section. 

The Poincare theorem is not applicable to the bunch of trajectories or even to two 
trajectories as it follows from the Lyapunov theory: if two points are situated in some small 
region of the phase space then they are not necessary come back to this region by moving 
along their trajectories. 

8.3 Mixing and weak limit 

The state p t = pt{x) on the compact phase space T is called mixing if its weak limit at 
t — > oo is a constant, 

lim pt(x) = const . 

More precisely, a dynamical system (T,tft, dp) has the mixing property [2H [32] if 

lim < /, U t g >= 

t^oo 

for every f,g e L 2 (T). Here U t g(x) = g((p t (x)). For the mixing systems the bunch of 
trajectories is spreading over the phase space, hence in the functional mechanics we have 
irreversibility. 

The method of the weak limit which generalizes the Poincare results and which can be 
applied to a wide class of dynamical systems is developed in [TTJ [12] . 

Connection with the irreversibility problem can be explained on the following example. 
Let us consider the function of two real variables 

F{t iP ) = e^f(p), 
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where f(p) is an integrable function. It is clear that the function F(t,p) is periodic in t if p 
is fixed and it has no limit as t — > oo. However, if we integrate the function F(t,p) over p, 

F(t) = J e^f(p)dp, 

then we get the function F(t) which already has the limit (by the Riemann - Lebesgue 
lemma) : 

lim F(t) = . 

t— >oo 



9 Dynamics of a Particle in a Box 

Dynamics of collisionless continuous medium in a box with reflecting walls is considered 
in [29], CHI EG] • This studied asymptotics of solutions of Liouville equation. In functional 
approach to mechanics, we interpret the solution of the Liouville equation as described the 
dynamics of a single particle. Here we consider this model in the classical and also in the 
quantum version for the special case of Gaussian initial data. 



9.1 Dynamics of a classical particle in a box 

Consider the motion of a free particle on the interval with the reflective ends. Using the 
method of reflections [52], the solution of the Liouville equation([H]) 

dp p dp 
dt m dq 

on the interval < q < 1 with the reflective ends we write as 



p(q, p,t)= V \p (q --t + 2n, p) + p (-q + —t + 2n, -p)] , (60) 
' m m 



where it is assumed that the function po has the Gaussian form 
One can show that for the distribution for coordinates 



Pc(q,t) = J p(q,p,t)dp (61) 

one gets the uniform limiting distribution (pointwise limit): 

lim p c (q,t) = 1 . 

t— >oo 

For the distribution of the absolute values of momenta (p > 0) 

Pa(p, t) = Pm(P, t) + Pm(-P, t) , 

where 

Pm(p,t) = / p(q,p,t)dq, 



as t — > oo we get the distribution of the Maxwell type (but not the Maxwell distribution): 

. . 1 . _(p~po) 2 (p+pq) 2 

hm p a {p,t) = —j=r\e + e P j . 
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9.2 Dynamics of a quantum particle in a box 

The Schrodinger equation for free quantum particle on the interval < x < 1 with reflecting 
ends has the form 

lf pA = _HL&± (62 ) 

dt 2m dx 2 

with the boundary conditions 

0(O,t) = O, <f){l,t) = 0, tel. 
Solution of this boundary problem can be written as follows: 

oo 

cP(x, t)= J*] [ip(x + 2n, t) - ip(-x + 2n, t)] , 

n=— oo 

where ip(x, t) is some solution of the Schrodinger equation. If we choose the function tp(x, t) in 
the form, corresponding to the distribution (|T9|) . then one can show that in the semiclassical 
limit for the probability density \(p(x, t)\ 2 the leading term is the classical distribution p c (x, t) 

dSU». 



10 Conclusions 

In this paper the functional formulation of classical mechanics is suggested which is based 
not on the notion of an individual trajectory of the particle but on the distribution function 
on the phase space. 

The fundamental equation of the microscopic dynamics in the proposed functional ap- 
proach is not the Newton equation but the Liouville equation for the distribution function 
of a single particle. Solutions of the Liouville equation have the property of derealization 
which accounts for irreversibility. It is shown that the Newton equation in this approach 
appears as an approximate equation describing the dynamics of the average values of the 
positions and momenta for not too long time intervals. Corrections to the Newton equation 
are computed. 

Interesting problems related with applications of the functional formulation of mechanics 
to statistical mechanics, to singularities in cosmology and black holes, and new interpretation 
of quantum mechanics we hope to consider in further works. 
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